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SYMBOLS  USED  AKD  UNITS  IN  EQUATIONS 


c  —  a  dimensionless  energy  ratio  of  matetlals 
e  --  unit  compressive  deformation,  ft.  per  ft. 

E  compressive  modulus  of  elasticity,  lb.  per  sq.  ft. 
f  —  stressing  turbulence  factor,  dimensionless 
H  —  length,  ft. 

S  --  unit  compressive  stress,  lb.  per  sq.  ft. 
t  --  time,  sec. 

u  —  internal  energy,  ft.  lb.  per  cu.ft. 

V  —  material  velocity,  ft.  per  sec. 

U  wave  velocity,  ft.  per  sec. 

X  --  distance,  ft. 

p  --  material  density,  slugs  per  cu.  ft. 


HYPER-ELASTIC  IMPACTS 


1.  THEORY  AND  ANALYSIS 

1.1  Introduction 

rHirtng  the  hlgh-veloclty  impacts  of  netal  objects,  stresses  beyond  the 
elastic  range  of  the  materials  are  developed  and  vaves  which  will  be 
designated  •'hyper -elastic  waves’*  are  propagated.  It  is  the  purpose  of 
this  report  to  investigate  such  one-dimensional  hyper-elastic  waves  in 
metal  rods  that  are  impacted  in  and  above  the  elastic  range. 

1.2  The  Nature  of  Hvper-Elastic  Waves 

It  seems  to  be  currently  popular  to  attempt  the  eKplanation  of  hyper¬ 
elastic  waves  by  contriving  an  equivalent  bulk  elasticity  fnr  metals  and 
treating  such-explosion  or  impact  produced  disturbances  as  classical 
fluid  waves.  Such  an  over-simplified  approach  fails  to  explain  some  of 
the  most  important  and  easily  observable  facts  of  such  a  disturbance. 

a.  It  is  not  a  single  wave  passing  through  the  material  since  the 
metal  near  the  impacted  end  of  a  bar  is  permanently  deformed  and 
"up-set"  when  the  metal  further  from  the  end  is  relatively 
unharmed . 

b.  A  strain  gage  record  of  the  disturbance  as  it  progresses  along  a 

_ ^bar  shows  it  not  to  be  a  classical  square  wave  front.  Instead, 

it  is  a  front  of  decreasing  steepness  with  definite  evidence  of 
inelastic  slip. 

c.  The  discontinuous  and  non-linear  nature  of  the  stress-strain  curve 
of  a  metal  does  not  suggest  that  its  load-deformation  nature  and 
action  can  be  described  by  a  single  number. 


It  therefore  seems  much  more  logical  to  base  the  analysis  of  such 
disturbances  on  the  observed  stress-strain  nature  of  actual  metals  rather 
than  to  treat  them  as  fluids  with  continuous  properties. 

1.3  The  Elastic  Recovery  of  Metals 

The  obvious  physical  properties  of  metals  are  summed  up  by  their  stress  S 
deformation  e  curves  which  characteristically  have  a  linear  portion  followed 
by  a  non-linear  section  terminating  at  the  ultimate  stress.  The  ultimate 
stress  in  compression  is  clearly  defined  for  brittle  materials  but  not  always 
so  evident  for  ductile  metals.  The  linear  portion  seems  not  to  be  altered 
by  the  rate  of  loading,  but  the  non-linear  portion  cm  be  changed  slightly 
by  the  rate  of  loading  during  testing,  particularly  at  hlgaer  strain*. 

The  non-linear  portion  is  referred  to  as  inelastir;  and  the  linear  portion 
as  the  elastic  range.  That  this  linear  range  is  not  truly  elastic  in  the 
reversible  sense  is  made  quite  evident  by  the  observed  fact  that  for  cyclic 
stresses  in  this  so-called  elastic  range  internal  hysteresis  losses  of 
considerable  magnitude  can  be  measured. 

Ano' her  indication  that  this  linear  and  so-called  elastic  portion  of 
the  stress-strain  curve  does  not  reveal  all  that  goes  on  occurs  when 
elastic  waves  arc  propagated  in  metal  bars.  For  a  metal  having  an  elastic 
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modulus  E,  ^  “  and  a  mass  density  p,  such  waves  should  propagate  at  a 
velocity .  However,  careful  tests  consistently  result  in  propagation 
velocities  greater  than  this. 

The  stress-strain  curve  depicted  by  Fig.  1,1  reflects  only  the  input 
work  required  to  deform  the  material  and  as  such  is  incapable  of  revealing 
what  became  of  this  input  deforming  work.  One  certain  way  of  determining 
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the  portion  of  this  Input  energy  that  Is  stored  elastically  Is  to  ask  the 
stressed  material  how  much  energy  it  can  return  when  the  stress  is  released. 
This  ' Jturnable  energy  is,  by  definition,  the  reversible,  elastically 
stored  energy. 

This  Is  accomplished  experimentally  by  stressing  the  metal  to  condition  A, 

Fig.  1.1,  and  then  determining  the  recovery  line  A-B  as  the  stress  is  released. 

If  this  is  done  at  various  stresses  up  to  the  ultimate,  each  time  on  a 

fresh  Specimen,  two  things  become  apparent.  The  recovery  lines  A-B  are 

straight  lines  and  all  of  the  recovery  triangles  A-B-C  are  similar  triangles. 

2 

eS  S 

If  In  each  case  the  recovery  energy  is  expressed  as  u  ■  -r”  or  ~  ,  the 

2  ‘'*“1 

recovery  modulus  Ej^  so  calculated  gives  the  recovery  curve  with  slope  Ej^ 
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as  Illustrated  by  Fig.  1.1.  Figures  1.2.  1.3.  and  1.4  show  such  test  results 
for  steel,  alumipum,  and  copper.  One  significance  of  the  recovery  modulus 
is  illustrated  by  Table  1.1. 


■■ 

■iHllli 

HH 

IHI' 

Steel 

15.15 

35.900.000 

29.300,000 

18,500 

16,700 

18,400 

Aluminum 

5.36 

12.400.000 

10,600,000 

18,200 

16,800 

17,900 

Copper 

17.28 

19.800.000 

16,000,000 

12,700 

11,540 

13,400 

TABLE  1.1 


As  is  shown  by  Table  1.1.  the  propagation  velocity  calculated  from  the  recovery 
modulus  Ej  comes  much  closer  to  checking  with  the  test  velocity  of  dis¬ 
turbances  than  does  the  value^E^/p  calculated  from  the  convent icnal  clastic 
modulus. 


1 . 4  Generalized  Stress-Strain  Curves 
Before  undertaking  the  analysis  of  hyper-elastic  waves,  it  is  helpful  to 
look  at  the  effect  of  the  stress-strain  curve  shape  on  the  propagation  of 
waves  in  materials.  It  was  first  shown  for  concrete  by  Talbot^  and  later 
for  other  materials  by  Gilkey  and  Murphy^  that  quite  different  appearing 
stress-strain  curves  can  be  reduced  to  essentially  similar  curves  for  each 
class  materials  by  making  a  unit  or  dimensionless  plot  of  S/S^  vs  e/e^ 
where  S  is  unit  stress  and  e  is  unit  strain.  Correspondingly,  and  e^ 


Talbot,  A.  N. ,  "Tests  of  Reinforced  Concrete  Beams,  Series  of  1905," 
Bui.  4,  Eng.  Ex.  Station,  University  of  Illinois,  Urbana,  1906. 

H.  J.  Gilkey  and  Glenn  Murphy,  "The  Percentage  Stress-Strain  Diagram  as 
an  Index  to  the  Comparative  Behavior  of  Materials  Under  Load,"  Iowa 
Engineering  Bulletin  159,  Iowa  State  College,  Ames,  1943. 
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UNIT  STRAIN  e  (10  in/in) 

Unit  stress  plotted  against  unit  strain  for  1018  steel  in  compression 
Specimens  of  3/4"  diameter  and  2”  long  were  used. 


UNIT  STRAIN  e  (1C  In/in) 

Unit  stress  plotted  against  unit  strain  for  24T4  alumlnuni  in  compression 
Specimens  of  3/4"  diameter  and  2"  long  vere  used. 
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UNIT  STRAIN  e(10”^  in/in) 

Fig.  1.4a.  Strain  energy  divided  by  the  unit  stress  times  unit  strain  (~^)  plotted 
against  unit  strain  (e)  for  copper  in  compression.  Specimens  of  3/4" 
diameter  and  2"  long  were  used. 


are  the  ultimate  values.  Such  curves  are  illustrated  by  Fig.  1*5. 


For  concrete,  Talbot  had  success  with  a  parabolic  equation  between 

S/S  and  e/e  but  to  cover  all  of  the  cases  of  Fig.  1.5,  a  more  versatile 
u  u 

equation  is  needed.  The  form  of  the  curves  of  Fig.  1.5  suggests  that  they 
should  be  described  in  three  variables.  Actually,  the  stress  S,  the  de¬ 
formation  e  and  the  total  stressing  energy  u  per  unit  volume  are  all 
Invr  ved  in  wave  action.  These  three  variables  can  reasonably  well  be 
related  by. 


u 

Se 


c 


(1.1) 


where  c  is  a  characteristic  constant  for  any  family  of  materials.  Since 
the  Internal  energy  is  changed  only  by  applied  stress, 

du  ■  Sde  (1.2) 
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.Equations  1  and  2  give 


u 

u 

u 


1 


u 


(1.3) 

(1.4) 

(1.5) 

(1.6) 


As  may  be  seen  from  Fig.  1.6,  £<>.  1.5  does  fit  steel  quite  well  in 

the  hyper>elastic  range.  Actually,  the  degree  to  which  Eq  1.5  can  made 

to  fit  test  stress-strain  curves  of  metals  is  largely  determined  by  the 

values  of  the  ultimate  stress  S  and  ultimate  deformation  e  that  are  used. 

u  u 

As  can  be  seen  from  Figs.  1.2,  1.3,  and  1.4,  the  ultimate  stress  in 
compression  is  not  clearly  defined.  However,  such  a  formulation  as  that 
of  Eq.  1.1  can  be  used  to  describe  both  linear  and  non-linear  materials. 


1.5  Waves  in  Linear  Materials 


f 


Fig,  1.6, 


Unit  stress  divided  by  ultimate  stress  plotted 
against  unit  strain  divided  by  ultimate  strain. 
Calculated  from  ^  C  and  dU  ■  Sde. 
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In  the  classical  case  of  a  stationary  unstressed  rod  struck  by  a 
sustained  impact  velocity  Vj^,  a  stress  wave  moves  through  the  material  at 
a  steady  velocity  Behind  this  front,  the  material  is  at  a  uniform 

stress  Sj  with  a  unit  deformation  e^  and  moving  at  Vj^.  The  wave  processes 
cu.ft.  per  sec.  of  material  and  the  transmitted  power  reaching  the 

front  must  give  this  material  its  internal  energy  u^  and  its  kinetic 

pvj 

energy  per  unit  volume 

r  “’'n 

L"i*— J 

The  stress  on  the  unit  area  must  cause  s luges  per  second  to 

gain  a  velocity  so 

Si  -  cMjVj  (1.8) 

The  velocity  Vj^  must  be  great  enough  to  take  up  the  slack  resulting  from  the 
unit  deformation  e^^  or 


Vi  -  ejWj 


For  a  linear  material  with  a  modulus  E, 


and 


Ii_ 

2E, 


(1.9) 


(1.10) 


(1.11) 
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From  Fig-  1.6  this  is  the  condition  for  a  material  with  a  linear  stress- 
strain  curve. 

Before  concluding  that  a  single  stress  wave  cannot  propagate  in  a  non¬ 
linear  material,  it  should  be  recalled  that  in  a  gas  with  the  non-linear 
curve  of  Fig.  1.5,  shock  waves  are  propagated.  It  is  also  known  that  the 
action  in  such  shock  fronts  is  highly  turbulent.  Turbulence  in  a  solid 
material  can  be  set  up  by  arguing  that  the  internal  energy  u  of  Eq.  1,1, 

Fig.  1.5,  and  Fig.  1.6  represents  only  the  final  energy  mechanically  stored. 
The  amount  of  supplied  energy  to  produce  this  would  be  greater  than  this 
fuj  where  f  is  a  number  greater  than  unity.  This  would  mean  that  energy 
[(l-f)uj  would  be  lost  as  heat  during  the  stressing  process. 

The  space  Eq.  1.9  would  not  be  altered 


V  >  e  U 

''l  I  I 


However,  the  energy  Eq.  1.7  would  be 


S^Vi  -Wj 


pv 

fui+ 


The  force  Eq.  1.8  would  remain 


and  the  non-linear  medium  would  be. 


(1.18) 


(1.19) 


(1.20) 


(1.21) 
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Combining  these  gives. 


WiPvJ  .WjpvJQe+l]  (1.22) 

Equation  1.22  can  be  satisfied  for  c  greater  than  1/2  only  by  a  f  less  than 
unity  and  this  cannot  be.  However,  for  c  less  than  1/2  (Fig.  1.6)  there 
is  a  turbulence  level  depicted  by  f  greater  than  unity  for  which  a  single 
wave  front  can  propagate.  Since  for  gases,  c  is  less  than  1/2  this  bears 
our  what  is  known  about  the  high  turbulence  in  gas  shock  waves. 

All  of  this  fits  the  following  physical  picture.  For  a  linear  material 
c  “  1/2  the  transmitted  energy  per  sec.  for  a  unit  area  which  reaches 

the  front  is  exactly  the  correct  amount  needed  to  supply  the  stress  energy 
and  the  kinetic  energy  involved.  Consequently,  such  a  wave  rc.eralule 
and  it  can  exist  either  as  a  compression  wave  or  as  a  recovery  or  rare* 
faction  wave.  For  a  material,  such  as  a  gas,  described  by  a  c  less  than 
1/2  the  energy  reaching  the  front  is  greater  than  that  needed  to 

supply  the  internal  energy  and  the  kinetic  energy.  Consequently,  if  there 
exists  a  mechanism  for  losing  energy,  (f  greater  than  unity)  a  compression 
wave  can  occur,  as  it  does  in  a  gas.  For  a  material  such  as  the  materials 
represented  by  c  greater  than  1/2,  the  amount  of  energy  transmitted  to 
the  front  Is  not  enough  to  supply  the  required  internal  energy  and  kinetic 
energy,  so  such  a  wave  cannot  occur. 

1 , 7  Hyper-Elastic  Waves 

For  any  analysis  of  hyper-elastic  impact  stresses  to  be  valid,  it 
must  be  capable  of  explaining  certain  physical  observations. 
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a.  When  strain  gag**  traces  are  determined  near  the  Impacted  end  of 


a  bar  A  and  at  the  center  B,  Fig.  1-8  Illustrates  typical  results. 
When  tfie  impacting  velocity  V  is  near  or  below  30  ft.  per  sec.,  two 
things  seem  apparent  on  such  films.  The  progress  of  the  "toe”  of 
the  wave  1  is  checked  by  the  elastic  velocity  using  the  recovery 
modulus  s  of  Table  1.1  and  the  full  deformation  point  2  is 

These  two  values 
of  E  are  explained  by  Fig.  l.l.  The  rise  time  1-2,  corrected 
for  strain  gage  length,  becomes  greater  as  the  strain  gage  is  moved 
down  the  bar.  This  suggests  that  such  an  elastic  wave  is,  in 
reality,  a  zone  and  not  a  front. 


checked  by  the  usual  elastic  modulus  as 
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b.  When  the  impact  velocity  V  is  raised  to  80  or  90  ft.  per  sec.,  the 
strain  gage  traces  of  Fig-  1.9  are  typical.  The  first  strain 
increase  1-2  remains  essentially  the  same  as  in  the  low  velocity 
case  of  Fig.  1.8.  However,  it  is  followed  by  a  later  very  large 
strain  which  quickly  goes  off  scale  for  gage  A  but  rises  more 
casually  at  position  B.  The  deformation  e^  is  consistently  at 
or  below  the  so-called  elastic  deformation  of  Figs.  1.2  to  1.4. 
The  rate  of  rise  of  the  following  hyper-elastic  deformation  2-3, 
Fig.  1.9,  becomes  progressively  less  along  the  bar. 


Time-microseconds 

Fig.  1.9 


The  hyper-elastic  wave  hypothesis  which  appears  to  be  consistent  with 
all  of  these  observations  is  illustrated  by  Figs.  1.10  and  1.11.  When  the 
bar  of  Fig.  1.11  is  Impacted,  an  elastic  wave  with  a  deformation  ej^ 
travels  along  the  bar.  This  is  followed  by  a  zone  in  which  the  elastic 
slip  e^  to  e^  occurs  and  the  back  of  this  zone  moves  down  the  bar  at  a 
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Fig.  1.10 


W  •  Wave  Velocity  V  »  Material  Velocity 
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Fig.  1.11 

velocity  W^.  Finally  this  is  followed  by  a  hyper-elastic  slip  from  e^  to 
e-j, which  brings  the  rod  velocity  up  to  the  Impact  value  V^.  At  lower 
impact  velocities,  this  last  action  can  be  absent. 
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1 . 8  Analysis  of  Hyper-Elastic  Impacts 
For  it  to  have  engineering  worth,  the  explanation  of  hyper-elastic 
impacts  outlined  in  Article  1.7  must  be  expressed  in  numbers.  Only  then 
can  it  be  verified  and  used. 

a.  Elastic  Wave  0- 1 .  Figures  1.10  and  l.ll 
From  Eqs.  1.10,  1.12,  and  1.13,  article  1.5 


(1.23) 

(1.24) 

(1.25) 


b.  Elastic  Slip  1-2.  Figures  1.10  and  1.11 

Test  results  illustrated  by  Figs.  1.8  and  1.9  suggest  that 
the  slip  deformation  1-2  actually  occurs  in  a  zone  at  finite  length. 
However,  to  keep  mathematical  complications  from  precluding  a 
workable  solution,  the  venerable  "lump  constant"  method  of  reducing 
it  to  an  equivalent  discontinuity  will  be  used. 

This  elastic  slip  wave  moving  at  a  velocity  overtakes 
(W^  -  Vj^)  cu.  ft.  per  sec.  of  metal  and  increases  its  velocity  from 
to  1^2*  A  force  balance  gives. 


(Force),  S2  -  -  p(W^  -  Vj)(V2  -  Vj) 


(1.26) 
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The  material  (W^  -  V^)  with  a  deformation  e^^  must  occupy  space 
^2  ”  '^2  *  deformation  or 

(Space)  (W^  -  Vj)(l  -  e^)  -  (M^  -  V^Xl  -  e^)  (1.27) 

The  transmitted  energy  per  sec.  reaching  the  slip  front 

must  supply  S^Vj  transmitted  to  the  elastic  wave  front  and  supply 
the  internal  energy  and  kinetic  energy  increase  for  (W^  -  Vj^)  cu. 
ft.  per  sec.  or 

r  p'^ii 

(Energy)  S2V2  -  SjVj  +  (W^  -  V^)  U,  -  u^  +  -5-  -  -f  ^  (1.28) 

While  the  slip  path  1-2  of  Fig.  1.10  is  not  acLt'ally  known,  a 
reasonable  approximation  is 

(S,  +  S  ) 

“2  ■  “1  “  ”4 - ^  ^®2  " 


and,  of  course, 


^2  "  ^2®2 


(1.30) 


Eliminating  from  Eqs.  1.26  and  1,27, 


(S-  -  S.)  V  (1-e  )-V  (1-e  ) 

«2  - ''l  * 


Similarly,  using  Eq.  1.26  to  eliminate  (W^  -  Vj^)  from  Eq.  1.28, 


(S,  -  S  ) 

^2'^2  -  ^1^  ^  wprri 


r 

-  vp  [“2  '  "'i 


2  2 

PV2  pVj- 

“T"  ■  '~2~ 


(1.32) 
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Using  Eqs.  1.23  and  1.30  to  eliminate  S^,  and  from  Eq.  1.31, 


„  ./IT.  .viiv-.Nr 

2  V  »  1  iPT  ('2  -  'l> 

»<^2  -  'ilf  r> 


—  (I  -  e,) 

0  2' 


(1.33) 


Using  Eqs.  1.23,  1.29,  and  1.30  to  eliminate  Vj^,  and  S2  from 


Eq.  1.32, 


h‘2^1  ■  ' 


riVili££ 

‘  ..  .X  L  2  <'2  '1) 


‘’<»2-'l\/p 


.d  fi!?] 

2  “  2  J 


(^.34) 


Clearing  Eq.  1.33 


■  ®1^^2  ’  ^  ^'^2®2  "  ®1®1^^®2  *  ®1^ 


!<l"®l)(V2"®l  ®iN/^  ^^"®2^^2"®1^^ 


)  (1.35) 


Collecting  terms  in  Eq.  1.35, 


p(l-e^)V2  -  V2  j^Sj^EjP  (l-e^)  +  p(l-  ®l/^  ^®2"®1^ 

+  e^E^Cl-e^)  -  (E2VVP^®2"®1^  ®1®1  ^®2  "  “  0  <1-36) 


v2  2\/S^e  V  I  ^1®1^^'®1^"^^2®2~^1®1^^®2'®1^ 
''2  ■  P  ^  2  p(l-e^) 


0  (1.37) 
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Clearing  Eq.  1.34 


2p(V^  -  )V2’'2  -  <"2  ■ 

+  (^2*^1*^  V2^^‘'2“®P  ^'^2 


(1.38) 


Collecting  terms  in  Eq.  1.38, 


pv^  [2V2  -  '2'2-'l'.]  f"l'2^2'[?«^'5^  ] 


2t[e\ 


-  (E2e2*^l*l) 


E'l'l  ^  V2><'2-"i>-'i'5] 


0  (1.39) 


or 


2EjeJ  - 

[^2®2'Ei®i] 

[<'2'2  ♦ 

E^.,)(p,-e,)-E,.2] 

P 

p2®2 

*  'I'l] 

-  0 

(1.40) 


Since  Eq.  1.37  combines  the  force  and  space  equations  and 
Eq.  1.40  combines  the  force  and  energy  equations,  their  simultaneous 
solution  describes  the  conditions  that  must  be  met  by  the  elastic 
slip  front,  1-2  Fig.  1.11,  to  enable  it  to  drive  the  elastic  wave  1 
ahead  of  it. 

Eliminating  V2  from  Eqs.  1.37  and  1.40, 


®1®1^^"®1^"^V2"^1®1^^®2"®1^ 

r  P(1  -  e^) 


2  3 

-^2®2"®1®1 

B®2'2  +  '1 

h>^®2-®l>  “  Vl] 

p 

Ve  *  h'l] 

(1.41) 


Clearing  Eq.  1.41 


EjBjd-ej)  1^262  +  EjcJ  -  < V2“®1®1^  ®1®1^^®2"®1^ 


(1.42) 

2Ejej(l-ej)-(E2e2-Ejej)(E2e2+Ejej)(e2-ej)(l-ej)+(E2e2-Eje^)E^eJ(l-e^) 


Simplifying  Eq.  1.42  gives 


(E2®2  ^1®1^^^2®2  “  ^1®1^^®2'®1^  " 


(E^e^  +  ‘  *^i®l^^®2  '  <1*^3) 

Equation  1.43  has  three  possible  solutions, 


(1  -  ep  -  1 

or 

®i-o 

(*2  -  Cj)  -  0 

or 

®2  "  ®1 

(1.45) 

Vz  ■  'I'l  ■  ° 

or 

S2  -  s, 

(1.46) 

Of  the  three,  e^  *  0  seems  to  fit  all  observed  conditions. 
When  e^  ■  0  is  used  in  Eq.  1.37,  the  result  is 


V 


2 


(1.47) 


Using  this  value  of  V2  and  =  0  in  Eq. 


This  value  of  in  conjunction  with  Wj^ 


1.31  gives, 


from  Eq. 


(1.48) 
1.24  goes 
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a  long  way  toward  explaining  the  fact  illustrated  in  Figs.,  1.8 
and  1.9  that  measured  deformations  bear  little  resemblance  to 
classical  square  wave  fronts. 

Remembering  that  the  "Iwap  constant"  method  here  used  is  really 
attempting  to  describe  what  is  probably  a  deformation  zone,  the 
following  picture  begins  to  make  sense.  If  the  front  of  this  zone 
travels  at  E^/o  and  the  rear  of  it  at  E^/p  .  the  zone 

will  get  progressively  longer  as  it  travels  along  the  bar  because 
Wj  is  greater  than  W^.  As  illustrated  by  Figs.  1.8  and  1.9,  this 
explains  why  the  deformation  front  becomes  longer  and  flatter  as 
the  disturbance  travels  along  the  bar. 

Hvper-Elastic  Slip  2-3 

The  hyper-elastic  slip  2-3  of  Fig.  I.IO  occurs  progressively 
over  a  varying  length  so  its  reduction  to  a  single  "lump  constant" 

effect  will  have  'to  be  based  on  average  conditions.  The  average 

e  e 
2  3 

deformation  in  the  slip  zone  may  be  taken  as  x — —  and  the 

V  +  V  ^ 

2  3 

average  material  velocity  in  the  zone  is  "■  .  The  average 

energy  increase  would  be  taken  as  half  of  the  maximum  value.  On 
this  basis  the  force,  space,  and  energy  equations  corresponding 
to  Eqs.  1.26,  1.27,  and  1.28  become 


Force 


Space 


Energy 


S3  -  -  P(W3  -  VjX-iy-J) 

S3V3  -  +  (H3  .  V3)  i  [P3.U3  +  ^ 


(1.49) 

(1.50) 

(1.51) 
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While  the  relationship  between  and  e^  is  determined  by  Eqs.  l.AO 
and  the  choice  of  one  of  Eqs.  1.44,  1.45,  or  1.46,  it  is  known 
that , 


so 

Force 


Space 


(W3-V2) (1-03/2-02,2)  -  CW3  -  V3)(l-e2) 


(1.54) 


2  ..2  „2 


1  r  ^^2 

Energy  S3V3  -  E2e2V2  +  (W3  -V2)  2  [“3  “  2  "  2  “  '  '5) 


Eliminating  from  Eqs.  1.53  and  1.54 


u  V  -  ''2<*-2/2  -  '’3/2> 

"3  ■  ’2  pCVj  -  Vj)  - 


(1.56) 


Using  Eq.  1.53  to  eliminate  (W3  -  V2)  from  Eq.  1.55, 


^3^3  “  ®2®2''2  p(V,  -  V, 


-"2°2>  r  ^ 

3  -  V2)  L“3  -  2  ■"  2  -  2  J 


(1.57) 


Clearing  Eq.  1.56 


PV2('^3"V2)(®3/2  ■  ®2/2^  ^  ^®3"^2®2^  ^®3/2"  ®2/2^ 


pV3(l-e2)(V3-V2)  (V3-V2) (1-62/2  ‘  ®3/2^ 
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Collecting  terms  In  Eq.  1.58 


(l-.j)cv^  .  (l.e^)2pv^v^  *  cv^(l-.j)  -  <S3-Ej.j)fe3-ej)  -  0  (1.59) 


2  PV2(l-e2^*^*3‘V2^^*1  "  *?> 

V  -  2V  V  +  — = _  i  £  Z  i  2 

3  ^  C(1  -  e,) 


(1.60) 


Clearing  Eq.  1.57 


'3^-''2>®3  ■  ''''2<»3-''2>V2  <*3-V2>t3  '  ^  ^ 


Collecting  terms  in  Eq.  1.61 


(S+  E2e2) 


3  2 


(1.61) 


"A  A 


-  i»3(S3«3.3)V3  +  pV^EjCj.  (Sj-E^Cj) (Uj-  -Sji  .  -2)  .  o 


or  22  (1.62) 

»  2  E  pV, 

2  2pV;E,e,  -  2(S,  -  E,e,)(u,-  -  — 2  ) 

3  ^  ^(SprEpp - -  (1.63) 


Eliminating  V  between  Eqs.  1.60  and  1.63 


2  2 

2  E  ej  pV^ 

7  P  a  _P  »  \ _  -±-A  _  _ 45 


PV2a-e)  -(S,-E2e^)(e,-e,)  2pV2E2e2-2  (S^-E.e^)  (u^-  -  -^) 

- - ^3^8202)  - ^  (1-6^) 


Clearing  Eq.  1.64 


pV2(l-e2)(S3.E2e2)  -  (S3-E2e2)  (63-62)  (S3  +  £262) 


(1.65) 
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or 


($^-£262)  (63-62)  (S3  +  £2*2^  “ 
0n6  solution  is 


S3  - 


£->6- 
2  2 


or 


V2 


(1.67) 


Which  is  the  case  of  no  hyper-elastic  slip.  This  leaves. 


—  —  '  (e3-e2)(S3  +  £262)  =  (1-02) (2U3  -  £262)  (1.68) 

Since  e2  will  not  exceed  .003,  the  (1  •  e^)  multiplying  factor 
may  safely  be  used  as  1.00  and  a  cubic  equation  avoided.  In  that 
way. 


(1.69) 


Using  V2  -  e^\f^ 
multiplying  factor 


from  Eq.  1.47  in  Eq.  1.60  and  again  letting  the 
(1-62)  be  replaced  by  1.00  gives. 


(1.70) 
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Similarly,  from  Eq.  1-56 


(1.71) 


(1.72) 


0.71) 


(1.74) 


(1.75) 


(1.76) 


(1.77) 

(1.78) 


1.9  Proof  of  Proposed  Test  Method 

The  analysis  of  Article  1.8  was  concerned  with  an  investigation  of  the 
elastic  and  sll?  waves  in  a  bar  struck  by  a  sustained  velocity  at  one  end. 
The  classical  simulation  of  this  In  the  elastic  range  is  to  have  a  bar  moving 
at  a  velocity  strike  an  Imagined  stationary  and  rigid  surface.  However, 
since  a  moving  bar  shot  from  a  gun  cannot  be  instrumented  with  strain  gages, 
another  approach  is  needed. 

The  ..method  used  was  to  strike  a  stationary  Instrumented  bar  squarely  with 
an  exact  image  bar  moving  at  a  velocity  of  2V^.  Before  this  approach  can  be 
accepted,  it  must  be  demonstrated  to  be  valid  for  the  case  of  the  elastic 
and  slip  waves  of  Article  1.8. 


Impactioi,  Bar 


wave  end 

Fig.  1.12 

This  will  be  proven  If  the  impacting  bar  of  Fig.  1.12  moving  at  has 
the  same  stress  waves  as  those  of  Article  1.8  when  the  velocity  of  the  impacting 
end  is  suddenly  reduced  to  -  V J2  by  impact  with  the  identical  instrumented 
bar. 

For  the  elastic  wave; 

Force:  S,  -  (W,  +  V  )p(V  -V,)  (1.79) 

— '  ■■  2  2  o  o  I 

Space:  (W^  +  V^Xl-e^)  *  (^2  ^^2^  (1.80) 
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Simplifying  Eq.  1.87  also  gives 


(1.90) 


(1.91) 


Now,  for  the  average  conditions  in  the  slip  zone  of  Fig.  1.12,  the  front 


moves  at  a  velocity  W.  and  the  hack  of  it  at  W  .  The  average  increase  in 
(V,  -  V  (e  +  e  ) 


velocity  is  ■■  '"  ^  and  the  average  unit  deformation  is 


Force ;  (S^-  S2)  “  c(Mj  +  V^) 


(V2  V3) 


(1.92) 


r  ®2m 

Space;  (W^  +  V^)  ll  -  -  "j-  ■  H  -  (W^  +  V^Xl  -  e^) 

1  r  '^'^3  '^'^2 1 

Energy:  S2V2  -  S3V3  -  (W3  +  V^)^  [“3  ‘  “2  ■*■  ~2 - T  J 


2  2 


(1.93) 


(1.94) 


Us  Ing 


®2  ■  '2*2 • 


*2*2  '2*2 


2  2 


(1.95) 


and  from  Eq.  1.90, 


-  •  '2'/? 


(1.96) 


In  Eqs.  1.92,  1.93,  and  1.94, 


Force:  S3-  £362 


(1.97) 


^3  ♦  V  'aV?]  E  -  •  <“3  *  ’'3>'>-'3> 

V2«'.-2l(?>  -  V3  *  I  [*3  ♦  V*2/?iE3  -  ^ 


(1.98) 


Eii^r^x:  E  e 


(1.99) 


Elimiiuitlng  from  Eqs.  1.97  and  1.98 


E  •  -  V3»-«2) 

^  ^•3-V 

I 

Uilng  the  value  of  -j  (Wj  + 


(1.100) 


'2*2  “*3 


-  f  ,  , . <»3-V2>L“3  •  -P  *  ¥  ■  -  V«2V^>!I 

2*2 '*0  M  *3*3  r  ir  T  W- 

»L*.  -  •  *3] 


The  deslfed  ^oof  can  be  established  by  solving  for 
Eqs.  1.  .'00  and  1.101. 

Doing  this  for  Eq.  1.100  gives. 


[''.-2'/?  -  *3] 


[*.  ■  •2'/?  -  'J'  • 


(e3  -  e^)(S3-E^ep 
0(1-62) 


(1.102) 


Similarly  from  Eq.  1.101, 


fv  ,/F,f  (*3-V2>«V*2°2> 

Lo  2Y  p  ^3]  p(S3  +  E,e  ) 


3  “2^2- 


(1.103) 
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EquatingEqsT  1.102  and  1.103  gives 


(63 -e^)  iS^-E^e^){2u^-E2ep 

Pd-e^)  “  ■ 


pCSj  +  E^e^) 


(1.104) 


Conceling  the  root  S3  -  which  is  that  of  a  wave  In  the  elastic  range 
and  not  of  present  concern,  gives 


2U3  - 
S3  +  E^e^ 


<*^3  "  *2> 

(I  -  *2^ 


(1.105) 


Since  this  is  exactly  the  conclusion  reached  in  Eq.  1.68  of  the  original 

^*3  ”  *2^ 

slip  wave  analysis,  the  value  of  ■  "  p(i-e  )  " 

used  in  Eq.  1.102  to  give 


V3  •  2V2V3  +  - 


[*,  -  -  h] 


(1.106) 


But,  from  Eq.  1.47 


-V. 


(1.107) 


So 


V3  -  e 


!\/?  ■  "o  -  '2n/? 


(1.108) 


or 


V  -  f 
3  2 


(iao9) 


which  Is  the  necessary  condition  for  the  validity  of  the  impacting  bar  scheme 
of  Fig.  1.12. 
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1.10  Interpretat ion  of  Analytical  Conclusions 

The  picture  that  has  been  contrived  is  that  of  the  front  of  an  Inelastic 
slip  zone  traveling  as  an  elastic  wave  at  velocity  E^/p  and  with  deforma¬ 
tion  e^'  imagined  rear  of  this  slip  zone  where  the  material  has  reached 

the  full  impacting  velocity  is  moving  at  a  velocity  W^,  which,  from  Eq.  1.76, 

Is  usually  much  less  than  W^.  This  slip  zone  of  varying  length  has  unit 
defortration  at  the  front ,  e^  at  the  back  and  was  calculated  by  assuming  a 
linear  average  of  the  kinetic  energy  and  Internal  energy  of  the  material 
in  the  slip  zone. 

As  was  described  by  Eq.  l.i»i»  and  illustrated  by  Fig.  1.10,  this  front 
traveling  elastic  wave  has  a  slightly  varying  thickness  as  it  travels  along 
a  rod.  The  front  of  it  is  traveling  at  the  rear  at 

where  Ej^  is  the  truly  elastic  or  recovery  modulus  of  tlie  v-alerlal  arl  E2  is 
the  commonly  used  el.'tstic  modules  of  the  material. 


Fig.  1,13 
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The  space  distribution  of  unit  defonnation  along  a  struck  rod  is  illustrated 
by  Fig.  1.13.  From  it,  the  time  variation  of  the  deformation  at  any  point  x 
may  be  deduced. 

The  first  disturbance  would  reach  x  in  sec.  and  by  sec.  the 

unit  deformation  would  be  e^-  At  any  greater  time  t  the  deformation  e  at 
any  x  Is 


or  in  ratios, 


p,*  ."I 

'e  'e  -1)1  -  - 

'32  ’Ll  L 


(c/c  )  *  1  . . 

'  '  2’x  <1  - 


(l.lll) 


X  2 

where  —  is  the  location  of  the  strain  gage  and  -- —  ,  the  location  of  e. 

,  la  2 

along  the  bar,  is  a  measure  of  elapsed  time. 


The  curves  of  Fig.  1.14,  calculated  from  Eq.  l.lll  for  strain  gages  at 
x/L  •  .05  and  x/L  ■»  .30  differ  from  observed  curves.  The  calculated  slip 
deformation  builds  up  more  rapidly  than  is  observed  in  tests.  This  indicates 


that  the  straight  line  deformation  distribution  of  Fig.  1.13  is  not  as  likely 
as  the  dotted  curve  indicated.  This  means  that  attributing  a  linear  average 
value  of  stress  energy  and  kinetic  energy  to  the  material  in  the  slip  zone 
gives  results  that  are  too  high.  This  would  imply  that  the  values  of  elastic 
deformation  62  calculated  from  such  a  "lump  constant"  analysis  should  turn 
out  to  be  a  bit  low  since  the  energy  left  behind  in  the  slip  zone  appears  to 
have  been  over  estimated. 
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LOCATION  OF  ELASTIC  FRONT 

*3 

Fig.  1.14.  Deformation  ratio  versus  location  of  elastic 

*2 

front  W^t/L.  Calculated  curves  are  compared 

with  the  observed  curves  for  the  following  assumed 
W  e 

....  1  .  .  3  slip  deformation  .  .. 

conditions:  rr-  ■  1.1,  — •  ■  "  *5.0 

W  ^2  ^2  deformation 

and  ~  ■  .2  for  one  gage  at  X/L  -  O.Ot  and  one 
**2 

at  X/L  -  0.3. 


1.11  Comparison  of  Analysts  and  T»8t8 

Since  the  elastic  deformation  e^.  Fig.  1.9,  can  be  determined  reasonably 

well  from  strain  gage  traces,  this  deformation  e^  Is  shown  for  steel,  aluminum, 

and  copper  in  Figs.  1.15,  1.17,  and  1.19.  The  theoretical  curves  were 

calculated  from  Eqs.  1.75,  1.74,  and  1.76  using  the  properties  of  Table  l.l 

^3  “3 

and  from  values  of  u,,  •=“  ,  and  taken  from  Figs.  1.2,  l.2a,  1.3,  1.3a, 

1.4,  and  1.4a. 

Values  of  Uj/S^e^  and  yield  values  of  e2  and  Eq.  1.74 

gives  corresponding  values  of  impact  velocity  needed  to  produce  such 
deformations.  The  velocity  of  the  back  of  the  slip  zone  was  calculated 
from  Eq.  1.76.  For  the  linear  range  of  the  stress-strain  curves  1.1,  1.3,  and 
1.4,  e^  “  «2’  slip  and  Th»B  linear  portion  of  the 

62  vs.  curve  is  shown  extended  on  Figs.  1.15,  1.17,  and  1.19.  Figures 
1.16,  1.18,  and  1.20  show  the  corresponding  values  of  62  and  plotted 
against  the  total  slip  deformation  e^. 

As  would  be  expected  from  Fig.  1.14,  it  was  at  times  difficult  to  know 
where  e^  ended  and  the  effects  of  the  slip  deformation  e^  began,  particularly 
for  strain  gage  locations  near  the  impact  end  of  the  bar.  As  was  shown  by 
Eq.  1.109,  the  velocity  of  Figs,  1.15,  1.17,  and  1.19  is  1/2  of  the 
velocity  of  the  shot  rod. 

Three  things  stand  out  in  this  comparison  of  test  results  with  analytical 
predictions: 

1.  Particularly  for  aluminum,  there  appears  to  be  distinct  evidence  of 
an  elastic  wave  with  no  apparent  slip  at  deformations  greater  than 

—^’thatof  the  elastic  limit  of  the  static  stress-strain  curve. 
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20  40  60  80  100  120  140  160  180 

VELOCITY  (ft/sec) 

1-15.  Strain  62  plotted  against  velocity  V3  comparing  the  theoretical 
curve  with  the  experimental  results  for  1018  steel. 
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STRAIN  (10  ’’  in/in) 

Pi-8*  1*18.  Curves  of  strain  03  against  strain  e_  and  the  wave  velocity  ratio  E2/P 

for  24T4  aluminum.  * 
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0  10  20  30  40  50  60  70 

VELOCITY  (ft /sec) 

Fig.  1.19.  Graph  of  strain  versus  velocity  V-  comparing 
the  experimental  data  with  the  theoretical 
curve  for  copper. 
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STRAIN  (10  in/in> 

^^8*  1.20.  Curves  of  strain  against  strain  e2  and  the  wave  velocity  ratio 
W^/fEj/p  for  copper. 


2.  While  there  is  considerable  experimental  uncertainty  in  the  experimental 
evaluation  of  the  deformation  e^,  the  consistency  of  the  trends 

'-'Supports  the  concept  of  elastic  waves  and  slip  zones  upon  which  the 
analysis  was  errected. 

3.  Figures  1.16,  1.18,  and  1.20  indicate  that  while  the  slip  deformation 
e^  increases  rapidly  with  the  striking  velocity  V^,  the  deformation 
02  of  fho  elastic  front  increases  but  slightly.  It  is  quite  signi¬ 
ficant  that  the  velocity  of  the  back  of  the  slip  zone  decreases 
rapidly  as  the  severity  of  the  impact  increases.  This  directly  affects 
the  ability  of  energy  to  escape  Into  the  material  at  the  point  of 
impact  and  thereby  increases  the  local  dar.iagc  rpsuitlng  fr"!::  the 
Impact . 
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II.  EQUIPMENT  AND  TESTING 


2.1  Static  Cowpresston  Tests 

Static  compression  tests  were  made  on  3/V  diameter  by  2"  long  specimens 
of  1018  steel,  24TA  aluminum,  and  copper.  The  tests  were  conducted  with  a 
60,000  pound  capacity  Olsen  testing  machine  and  a  microformer  strain-load 
recorder  (a  drum  stylist  type  recorder).  A  specimen  was  loaded  to  a  certain 
stress  and  then  relieved  with  the  load  and  strain  being  recorded  during  both 
the  stresslngand  relieving  periods.  A  new  specimen  was  used  for  each 
different  load  to  eliminate  work  hardening  effects.  Strain  gages  were  also 
used  simultaneously  to  check  the  strain  an.!  stresses  in  the  elastic  region 
with  those  indicated  by  the  recorder.  These  tests  were  run  because  ci  the 
lack  of  suitable  compression  data  for  these  matcrlula  and  also  because  the 
dynamic  tests  were  to  be  performed  under  fast  compression  loading.  These 
data  were  used  in  the  development  of  the  theory  of  the  dynamic  part  of  the 
experiment.  The  actual  stress  strain  curves  of  1018  steel,  24T4  aluminum, 
and  coppe^ar^  shown  in  Figs.- 1.2,  1.3,  and  1.4  respectively. 

It  1 ^  definitely  seen  from  the  recorded  information  that  the  recovered 
energy  for  any  applied  load  is  always  less  than  the  input  energy.  This  is 
true  even  in  the  elastic  region.  Thus  some  permanent  strain  is  always 
present  and  the  slope  of  the  recovery  line  is  always  steeper  than  the 
previous  loading  line.  Stated  otherwise,  the  recovery  modulus  of  elasticity 
It  always  greater  than  the  loading  modulus  of  elasticity.  Figure  1.1  illustrates 
this  phenomena. 
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2.2  Dynamic  Test  Apparatus 


A  8p«jclnl  gun  consisting  of  one  solid  section,  a  breech  solenoid  section, 
and  a  barrel  extension  alignment  section  vas  designed,  fabricated,  and  used 
to  accelerate  the  metal  rods.  It  has  a  smooth  bore  of  O.SOO  inches  in 
diameter  and  Is  chambered  to  fire  a  50  caliber  cartridge.  The  gun  is 
mounted  on  two  steel  blocks  whose  design  allows  motion  in  the  axial  as 
well  as  the  transverse  direction  for  alignment  purposes.  The  mountings  are 
bolted  in  place  to  a  stationary  metal  table  which  holds  the  entire  test  set 
up,  and  the  table  is  rigidly  fastened  to  the  wall  of  the  concrete  tunnel  in 
which  the  tests  have  been  conducted.  An  air  gun  was  also  developed  by 
adapting  a  gas  solenoid  valve  to  the  already  existing  gun  barrel  and  it  is 
used  for  low  rod  velocities.  Photographs  of  the  powder  gun  and  air  gun  are 
shown  in  Figs.  2.1  and  2.2  respectively.  A  schematic  diagram  of  the 
complete  test  set  up  is  shown  in  Fig.  2.3. 

2.3  Test  Rods 

Twelve  inch  long  rods  of  1018  steel,  24T4  aluminum,  and  copper  were 
machined  to  as  close  to  0.499  inches  in  diameter  as  possible.  The  ends  were 
machined  and  all  machine  marks  were  removed  with  an  aluminum  oxide  cloth.  This 
was  done  in  order  to  reduce  to  a  minimum  any  concentrated  local  stresses 
which  might  be  set  up  upon  Impact.  Of  every  pair  of  rods  machine,  one  rod 
was  to  be  accelerated  in  the  previously  mentioned  50  caliber  gun  while  the 
other  was  to  be  used  as  a  stationary  rod  instrumented  with  strain  gages  at 
various  locations.  The  back  end  of  each  stationary  rod  was  placed  firmly 
against  a  hardened  and  rigidly  fastened  4340  steel  surface  so  that  little 
or  no  deformation  of  the  rod  at  that  location  would  exist,  thus  creating  a 


~  48 


Fig.  2.3.  Schematic  Diagram  of  Test  Setup. 


fixed  end  condition.  The  dimension  of  12  inches  for  each  rod  Insured  length 
enough  to  measure  the  propagated  velocity  of  the  wave  front  with  considerable 
precision.  Also,  because  about  5  inches  of  the  moving  rod  was  still  in  the 
gun  barrel  upon  impact,  good  alignment  was  maintained.  The  impact  and 
stationary  rods  were  aligned  before  each  firing  so  as  to  impact  squarely  and 
uni formly. 

The  12  inch  long  metal  rods  were  first  fired  into  a  box  of  sand  without 
impacting  with  other  rods  so  as  to  determine  either  the  mass  of  powder  or 
air  pressure  (depending  upon  which  gun  was  used)  required  to  produce  the 
subsequent  desired  velocities.  The  velocity  measuring  system  developed  and 
utilized  is  discussed  in  a  later  section  of  this  report.  Velocities  as  high 
as  1700  ff^/sec  were  measured,  but  it  was  later  found  out  that  £rc^  10  ft /sec 
to  500  ft /sec  was  the  most  practical  range  in  whicr.  to  vor<<  and  record  good 
data.  Graphs  of  powder  load  against  rod  velocity  and  air  pressure  against 
rod  velocity  are  shown  for  steel  only  in  Figs.  2.4  and  2.5  respectively. 

These  curves  were  used  only  to  predict  the  mass  of  gun  powder  or  air  pressure 
needed  to  produce  a  desired  velocity  for  the  velocity  was  measured  on  each 
firing. 

2.4  Dynamic  Strain  Measurements 

Budd  and  SR-4  strain  gages,  designed  especially  for  dynamic  loading 
conditions,  were  placed  at  various  locations  on  the  stationary  rod.  Each 
strain  gage  was  wired  into  a  separate  amplifier  designed  specifically  for 
the  strain  gages  and  then  connected  directly  to  one  or  more  oscilloscopes. 

The  amplifiers  and  strain  gages  are  discussed  In  Section  3  of  this  report. 

The  vertical  scale  of  the  oscilloscope  in  volts  per  cuntlmeter  was  converted 
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Fig.  2.4.  Graph  displaying  veloc; 
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to  inches  per  inch  strain  per  every  cenLimeter  through  an  equation  developed 
for  each  amplifier  and  which  is  also  discussed  in  Section  3.  Therefore,  the 
photographs  taken  by  the  cameras  attached  to  the  oscilloscopes  were  essentially 
of  strain  versus  time.  Sketches  of  some  of  the  types  of  photographs  taken 
are  illustrated  in  Figs.  2.6  and  2.8.  The  strain  labeled  e^  in  a  previous 
discussion  is  actually  the  vertical  deflection  of  the  almost  vertical  part 
of  the  trace  shown  in  Fig.  2.8a.  A  different  e^  was  indicated  and  measured 
for  each  different  impact  velocity.  This  experimental  data  of  one  half  the  rod 
Impact  velocity  versus  strain  can  be  seen  in  Figs.  1.15,  1.17,  and  1.19 
for  1018  steel,  24T4  aluminum,  and  copper  respectively,  and  there  compared 
with  the  theoretical  curves.  The  actual  data  from  which  the  experimental 
curves  were  constructed  appear  in  tables  A.  1,  A. 2,  and  <  of  the  Aj  pendiv 
for  steel,  aluminum,  and  copper  respectively. 

The  wave  velocity  was  measured  by  placing  two  strain  gages  on  a 
rod  a  considerable  known  distance  apart  in  the  same  plane  of  reference 
and  observing  the  time  in  micro  seconds  between  which  the  front  and  back 
strain  gages  were  first  disturbed.  The  average  of  the  wave  velocities 
for  all  three  materials  are  reported  in  Table  1.1.  The  wave  velocity  Wj^ 
was  found  experimentally  to  be  Independent  of  how  hard  the  bar  was  struck 
and  thus,  for  any  particular  material,  constant  for  all  Impact  velocities. 

One  interesting  observation  made  which  further  strengthens  the  theory 
previously  discussed  is  the  fact  that  a  strain  very  close  to  the  Impact  end 
displayed  almost  a  vertical  elastic  front,  whereas  the  slope  continued  to 
decrease  for  strain  gages  placed  farther  away  from  the  Impact  end.  This 
indicates  strongly,  as  mentioned  in  Section  1.7  that  because  travels 
faster  than  ^2  distance  between  them  becomes  greater  as  they  move 
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dovm  the  rod,  and  the  rise  tinie  for  the  elastic  wave  front  or  rone  Is 
Increased.  This  is  Illustrated  in  Fig.  2.6. 


Fig.  2.6.  Sketch  displaying  rise  time  of  wave 
front  versus  strain  gage  position. 

Because  and  arc  constant  values  for  a  particular  material,  theory 
predicts  the  rise  time  of  the  elastic  slip  zone  to  be  linearly  increasing 
with  distance.  A  graph  displaying  the  theoretical  rise  time  versus 
distance  is  shown  in  Fig.  2.7  and  there  compared  with  some  experimental 
values  for  1018  steel  only. 

One  other  item  that  should  be  mentioned  concerns  that  of  recognizing 
the  correct  strain  gage  output  voltage  which,  after  being  amplified,  was 
photographed  as  it  appeared  on  an  oscilloscope.  This  voltage  value  was 
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used  directly  in  the  calculation  of  strain  as  is  discussed  in  Section  3  of 
this  report.  Three  different  types  of  photographs  of  the  accumulated  data 
which  required  careful  Judgment  to  decipher  arc  shown  in  Fig.  2.8.  One 
type  that  was  not  so  difficult  to  interpret  is  that  of  Fig.  2.6. 


Fig.  2.8.  Sketch  displaying  different  type  photographs 
taken  of  the  strain  gage  output. 

For  the  type  of  photograph  Illustrated  in  Fig.  2.8a,  the  distance  A  was 
taken  to  be  the  correct  voltage  value  because  it  represented  an  average 
strain  of  the  material  processed  by  the  wave  front  W^.  For  the  type 
of  photograph  Illustrated  in  Figs.  2.8b  and  2.8c  .distances  B  and  C 
were  taken  to  be  the  correct  voltages  simply  because  they  represent  the 
linear  range  of  the  strains.  Other  types  of  photographs  taken  presented 
little  or  no  difficulty  to  decipher. 

Finally,  of  the  two  kinds  of  gages  used,  the  SR-A  strain  gage  generally 
Indicated  strain  values  which  fell  very  close  to  the  theoretical  curves  as 
compared  to  the  Budd  strain  gages  whose  strain  Indications  fell  further 
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away.  So  far,  no  explanation  for  this  difference  has  been  found,  and  all 
data  taken  from  both  gages  arc  displayed.  The  majority  of  the  SR-4  gages 
were  used  on  steel  and  very  few  on  copper  and  aluminum;  therefore,  the 
experimental  points  which  appear  very  near  the  theoretical  curves  are  more 
abundant  on  the  steel  curve  than  on  the  others. 


3 .  INSTRDMESIAT lOK 


3.1  Defornation  Measureiaents 
3.1.1  Strain  Transducers 

Recent  literature  has  shown  that  a  number  of  methods  have  been  used  to 

1,2,3  U 

measure  deformations  in  materials.  •  of  these  methods  the 

strain-gage  technique  appears  to  bo  the  simplest,  cheapest,  and  the  least 
time  consuming  and  was  therefore  adopted  for  this  project.  The  strain 
gage  instrumentation  is  somewliat  sirplified  in  this  case  since  most  of  the 
phenomena  to  be  measured  occur  within  200  microseconds  after  the  rod  impact. 
For  such  an  application,  it  is  not  necessary  to  use  temperature  compensated 
strain  gages  or  dummy  gages,  a  single  gage  attached  to  the  rod  at  the 
desired  spot  and  in  the  desired  orientation  being  sufficient  to  prov^  Je  strain 
measurements. 

Figure  3.1  shows  the  strain  gage  circuit  used  with  both  the  SR-4  and 

♦ 

the  Budd  type  strain  gages.  The  characteristics  of  these  gages  are  also 
given  in  the  figure.  Both  the  SR-4  and  the  Budd  gages  are  of  the  iso¬ 
elastic  type,  the  former  being  a  wire  gage,  the  latter  an  etched  foil  gage. 

^  Walsh,  J.  M. ,  Rice,  M.  H. ,  McQueen,  R.  C, ,  and  Yarger,  F.  L. ,  Physical 
Review.  Vol.  108,  p.  196,  Oct.  1957. 

2 

Bell,  J.  F. ,  Journal  of  Applied  Physics,  Vol.  27,  No.  10,  p.  1109,  Oct.  1956. 

3 

Curtis ,  C.  W. ,  Int .  Symp.  on  Stress  Wave  Propagation  in  Materials .  Ed.  by 
N.  Davids,  Interscience  Publishers,  New  York,  1960,  p.  15. 

4 

Elchelberger ,  R.  J. ,  Int .  Symp.  on  Stress  Wave  Propagation  in  Materials .  Ed. 
by  N.  Davids,  Interscience  Publishers,  New  York,  1960,  p.  15. 

* 

Strain  gages  used  were  Budd,  type  321-160,  obtained  from  Budd  Company, 

P.  0.  Box  245,  Phoenlxvllle ,  Pa.,  and  SR-4  type  C-11  obtained  from 
Baldwin-Llma-Hamilton  Corp. ,  Electronics  and  Instrumentation  Div.,  Waltham, 
Mass. 
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SR-4 


Budd 


Type 

R 

g 

Gage  Factor 


Active  Gage 
size 


C-11 
300  t  2 
2.86  t  2% 
100  +  I 
12 

1/8  X  1/8 


Fig.  3.1.  Strain-Cage  Circuit  for  Use  with 
SR-4  and  Budd  Cages. 


321-160 
160  -f  .2  ohms 
2.98  +  \X 
100  1  ohms 

6  volts  (approx.) 
1/8  X  1/8  inches 


These  gages  were  chosen  because  of  their  small  size,  high  gage  factor,  and 
availability.  SR-4  gages  can  be  fastened  in  place  by  simply  gluing  the  gage 
to  the  test  rod  with  either  Duco  or  SR-4  solvent-type  cernent.  Theii*'  cements 
require  approximately  12  hours  drying  time  before  use.  The  Euud  gages  must 
be  fastened  on  with  an  epoxy-type  cement  such  as  Budd  GAS  epoxy  adhesive. 
This  epoxy  needs  an  elevated  temperature  for  proper  curing,  a  typical  curing 
rate  being  200®  F.  for  one  hour. 

The  output  of  the  strain-gage  circuit  of  Fig.  3.1  can  be  calculated  as 
a  voltage  which  is  directly  proportional  to  the  strain  on  the  gage.  The 
output  voltage  E  is 


E 

o 


(3.1) 


If  R 

g 


changes  by  ^ 


then  changes  by  AE^  so  the  output  is 


E  +  AE  - 
o  o 


E 

b  Vr,  +  ..  iR^; 


(3.2) 


61 


3.2  Amplification,  Matching,  and  Recording 

A  block  diagram  of  the  firing  tunnels  and  control  consol  is  shown  in 
Fig.  3.2.  As  shown,  the  output  of  each  strain  gage  circuit  is  carried  to 
an  oscilloscope  through  an  amplifier,  a  cathode-follower,  and  a  terminated 
co-axlal  cable.  Channels  4,  5,  6,  and  7  are  Identical  in  construction  and 
nearly  identical  electrically,  differing  only  slightly  in  their  gains.  The 
overall  gain  of  each  channel  is  tabulated  in  Fig.  3.2.  This  gain  represents 
the  voltage  gain  between  the  input  to  the  amplifier  and  the  output  at  the 
end  of  the  terminated  co-axial  cable.  The  complete  circuit  diagram  for  one 
of  these  channels  is  shown  in  Fig.  3.3. 

The  velocity  with  which  an  elastic  wave  travels  down  a  steel  rod  is 
approximately- 17 ,000  feet  per  second.  Assuming  the:  •  I’l*?!!'-  Wei  -p  to  have  a 
very  sharp  front,  it  will  then  take  approximately  .6  inlf'rnspcnnti’’  fci  this 
wave  front  to  pass  beneath  a  strain  gage  whose  axial  length  is  1/8  of  an  inch. 
For  this  reason  the  strain  gage  amplifiers  are  designed  to  have  a  rise  time 
of  0.1  microsecond.  Figure  3.4a  Is  a  photograph  showing  a  rise  time  of  .3  ^secs 
which  is  typical  of  channels  4,  5,  6,  and  7,  This  response  time  is  that  of 
the  complete  system  including  the  strain  gage  amplifier,  cathode-follower, 
co-axlax  cable,  termination,  and  oscilloscope. 

The  output  of  the  strain  gages  during  a  test  Is  photographed  directly 
from  the  oscilloscope  using  Polaroid  Land  cameras.  A  typical  photograph  of 
strain  gage  outputs  in  shown  in  Fig.  3.4b.  All  strain  gage  channels  may  be 
used  simultaneously.  For  purposes  of  identification,  the  strain  gage  amplifier 
channels  are  numbered  channels  4,  5,  6,  and  7. 

3 . 3  Velocity  Measurements 

The  velocity  of  the  rod  as  It  approaches  Impact  with  the  fixed  rod  must 
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Fig.  3.3.  CIRCUrr  DIAGRAM  FOR  CHANNEL  4.  CIl^JJNELS  5,  6,  AND  7  HAVE  THE  SAME 
DIA(3tAM  BUT  DIFFER  SLIGHTLY  IN  PHYSICAL  ARRANGEMENT  AND  GAIN 


be  measured.  In  many  experiments,  the  velocity  o£  the  projectile  can  be 
measured  by  allowing  it  to  pierce  two  separate  grids,  spaced  a  fixed  distance 
apart.  The  time  necessary  for  the  projectile  to  travel  from  one  grid  to  the 
other  can  be  precisely  measured  and  the  velocity  calculated.  Such  a  method, 
however,  is  not  desirable  when  a  precision  impact  is  necessary,  especially 
when  using  rods  with  accurately  machined  ends.  It  is  possible  that  fragments 
of  the  grid  may  be  carried  ahead  of  the  traveling  rod.  Such  fragmentary 
pieces,  depending  on  the  material  used  for  the  gird,  may  cause  distortions 
In  the  Impact . 

The  velocity  measuring  device  to  be  described  here  is  an  optical  device 
eliminating  the  possibility  of  any  grid  fragments  distorting  the  Impact  and 
also  removing  the  possibility  of  any  external  applied  forces  on  the  rod  as 
It  moves  toward  the  Impact.  Figure  3.5  is  a  simplified  urawing  cT  this  system. 
Light  from  the  lamp  is  collected  by  lens  A  and  sent  parallel  through  the 
elongated  silt.  The  fixed  rod  is  positioned  such  that  the  impact  end  of  it 
Interrupts  a  slight  amount  of  the  light  at  the  impact  end  of  the  slit.  The 
narrow  beam  of  light  is  then  collected  by  lens  B  and  focused  on  a  photo  cell. 
The  output  of  the  photo  cell  is  fed  through  a  cathode  follower  and  along  a 
matc.ied  cable  directly  to  an  oscilloscope.  As  the  traveling  rod  leaves  the 
gun  and  moves  toward  the  impact  end  of  the  fixed  rod,  It  Interrupts  more  and 
more  of  the  narrow  beam  of  light.  The  output  of  the  photo  cell  begins  to  drop 
as  soon  as  the  traveling  rod  enters  the  light  beam  and  continues  to  drop 
linearly  as  the  traveling  rod  approaches  the  impact  end  of  the  fixed  rod.  After 
Impact  occurs,  there  is  no  further  change  in  the  photo  cell  output  as  long  as 
the  two  rods  stay  together.  If  after  some  length  of  time  the  traveling  rod 
bounces  back  out  again,  this  will  be  indicated  by  an  Increasing  voltage  at 
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Fig.  3.5.  SIMPLIFIED  DIAGRAM  OF  OPTICAL  VELOCITY  MEASURING  DEVICE 


the  photo  cell  output.  Ftora  a  photograph  of  the  photo  cell  output  as  recorded 
on  an  oscilloscope,  the  time  necessary  for  the  moving  rod  to  interrupt  the 
light  beam  can  be  measured.  This  tine,  and  the  measured  vldth  of  the  beam 
being  interrupted  (d  in  Fig.  3.5),  gives  the  average  velocity  of  the  rod 
just  prior  to  impact.  It  is  also  possible  to  determine  how  long  the  two 
bar direct  contact  with  each  other  and  the  nature  of  any  bounce  that 
may  take  place.  A  more  detailed  drawing  of  the  velocity  box  is  shown  in 
Fig.  3.6. 

The  cathode  follower  associated  with  the  photo  cell  is  mounted  inside 
the  box  immediately  beside  the  photo  cell.  The  photo  cell  circuit  is  sViown 
in  Fig.  3.7.  A  standard  12  volt  automobile  battery  is  used  to  power  the  12 
volt  10  watt  automobile  globe  used  as  a  light  source.  Halt- Incl.  plexi¬ 
glass  windows  are  provided  on  either  side  of  the  rnpncl  area  to  protect  the 
lenses.  The  unit  is  placed  near  enough  the  end  of  the  gun  such  that  impact 
Is  achieved  before  the  traveling  rod  completely  leaves  the  gun.  In  this 
way,  the  moving  rod  is  held  firmly  in  place  and  a  more  accurate  impact  can  be 
obtained.  A  typical  velocity  picture  for  a  steel-to-steel  impact  at  medium 
velocities  is  shown  in  Fig.  3.8. 
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Fig.  3,7.  PHOTO  CELL  AND  CATHODE  FOLLOWER  CIRCUIT  DUGRAM 
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Fig.  3.8.  Velocity  trace  for  steel-to-steel  impact 
showing  some  bounce.  Sweep  .5  m  sec/cm., 
amplitude  .2  v/cm,  d  =•  2.21",  t  ■  .85  m  sec. 
velocity  =  217  ft/sec. 


3.4  Tr iRRcr  C i rcii its 
3,41  Deformation  Trigger 

The  trigger  signal  used  to  start  tlie  sweep  circuit  on  those  oscilloscopes 
displaying  the  deformation  voltages  is  generated  at  the  time  of  Impact.  The 
method  used  to  obtain  this  trigger  signal  can  best  be  explained  by  referring 
to  Fig.'S.^.  As  shown,  the  gun  barrel  is  completely  insulated  from  ground 


> 

> 

> 

>  loon 


To  scope 


I 


.01  gfd 


Fixed  Rod 


Back-up 
✓  Block 


Fig.  3.9.  Deformation  trigger  circuit. 
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and  Is  initially  at  a  potential  of  approximately  45  volts  as  determined  by 
the  resistor  voltage  divider.  The  O.Ol  microfarad  capacitor  is  also  charged 
to  45  volts.  The  fixed  rod  is  grounded  as  near  to  the  impact  end  as  possible. 

At  the  time  of  impact,  the  traveling  rod  is  still  in  the  gun  and  therefore 
when  It  hits  the  fixed  rod,  the  gun  is  grounded.  The  capacitor  then  discharges 
with  a  theoretical  time  constant  of  1  microsecond.  However,  due  to  small 
inherent  inductances  in  the  system,  some  overshoot  is  present.  The  oscilloscopes 
are  adjusted  to  trigger  as  soon  as  the  voltage  drops  by  approximately  20  volts. 
Triggering  Is  reliable  within  approximately  one  half  of  a  microsecond  and  as 
many  as  five  oscilloscopes  can  be  triggered  simultaneously  with  this  system. 

3.42  The  Velocity  Trigger 

A  separate  trigger  circuit  is  necessary  to  p.'ovide  a  Signal  to  start 
the  sweep  on  the  velocity  oscilloscope  at  some  time  ju.st  prevlo 's  to  the 
instant  the  traveling  rod  enters  the  light  beam.  This  is  accomplished  by 
stretching  a  very  fine  cotton  thread  across  the  hole  in  the  velocity 
measuring  box  through  which  the  traveling  rod  enters.  This  cotton  string 
Is  in  turn  attached  to  a  number  34  copper  wire  which  forms  a  part  of  the 
trigger  circuit  as  shown  in  Fig.  3.10. 


Fig.  3.10.  Velocity  Trigger  Circuit 


This  small  copper  wire  forms  a  short  across  one  of  the  27  K  resistors  forcing 
the  trigger  output  to  be  at  0  potentials.  As  the  traveling  rod  enters  the 
hole,  it  strikes  the  cotton  cord  which  in  turn  breaks  the  number  3^  wire 
strap.  The  cutput  then  immediately  arises  to  approximately  -fAS  volts,  and 
the  scope  triggers.  The  cotton  string  is  adjusted  such  that  it  breaks  at 
the  same  tlmcthe  copper  wire  breaks  and  no  fragments  are  carried  into  the 
impact  area. 

In  cases  where  only  velocity  measurements  are  being  made,  it  is  possible 
to  place  the  copper  wire  strap  directly  over  the  hole  itself  eliminating  the 
need  for  the  cotton  string.  This  cannot  be  done  when  deformation  measurements 
are  being  made,  however,  because  the  copper  strap  is  at  ground  potential  and 
will  give  a  premature  deformation  trigger  "ignal  at  rhe  instant  the  traveling 
rod  hits  and  breaks  the  copper  wire  strap.  For  conv^niehci.-,  the  trtj.ter- 
circuit  batteries  and  components  are  mounted  insldn  veloelty  measuring 
box  in  such  a  way  as  not  to  interfere  with  the  light  beam. 


3,5  Accuracy  of  the  Battery  and  Strain  Gage  Circuit 
The  gage  resistance  of  a  Budd  strain  gage  is  160  +  .20.  This  gives 
an  accuracy  figure  of  +  0.1257..  The  ballast  resistor  of  the  strain  gage 
circuit  (See  Fig.  3.1)  is  a  17.  resistor.  The  battery  voltage,  being  measured 
with  a  calibrated  instrument,  can  be  read  with  an  accuracy  of  -t  17..  The 
change  in  the  output  voltage  of  the  strain  gage  circuit  of  Fig.  3.1  is  given 
by  Eq.  3.4  which  la  repeated  here. 


(3.4) 


where  a  has  been  defined  as  R. /R  . 

b  g 


For  the  errors  of  Rj^  and  R  as  given  above, 
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the  m.:ixinium  error  in  a  will  be  +  1.125”.  The  tern  /R  has  a  maximuin 

eg 

percentage  error  of  0,257„.  Of  all  of  the  terns  in  the  tierK'vninat or  of  the 

quantity  contained  within  the  second  set  of  brackets  in  Eq.  3.4,  the  term  a 

has  the  largest  percentage  error.  Thus,  the  largest  maximum  error  that  can 

occur  in  the  sum  of  terms  1  +  a  +  (ZR  /R  )  is  1.1257..  Equation  3.4  can  be 

g  g 

rewritten  indicating  the  percentage  error  in  each  of  the  terms  as  shown  in 
Eq.  3.10. 


AE 


E,(t  n)(YT 


-)(+  2.257.) 


/R 

R 


(*  0.257.) 


(f  1.1257.) 


Ll  +  a  + 


(3 


This  makes  the  net  maximum  error  in  the  outnut  voltage  of  the  strain  gage 

circuit  AE  equal  to  about  +  4.637.. 
o 

3.51  AmpH f ier  and  Osc i 1 loscopc  Cal ihrat ion  Accuracy 
The  gain  of  each  amplifier  was  measured  from  the  input  to  the  amplifier 
to  the  final  scope  termination  and  also  included  the  oscilloscope  amplifiers 
themselves.  An  input  wave  was  placed  directly  on  the  input  of  the  amplifier 
and  also  compared  directly  to  the  amplified  wave  on  the  oscilloscope.  It 
is  estimat  U  that  the  accuracy  of  this  overall  calibration  is  of  the  order 
of  -t  37o.  This  includes  the  amplifier,  cable,  cable  termination,  oscilloscope 
amplifier  and  reading  of  the  oscilloscope  face.  Since  the  output  voltage 
of  the  strain  gage  circuit  must  be  multiplied  by  the  amplifier  gain  to  obtain 
the  total  output,  the  errors  of  tlie  two  quantities  must  bo  added  giving  a 
total  overall  maximum  deviation  of  +  7.257.  It  should  be  pointed  out,  however, 
that  in  each  case  we  have  taken  Llic  worst  possible  error.  It  is  possible 


.10) 
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that  consistency  In  leaflinp,  data  and  ncasiiring  distances  on  the  oscilloscope 
face  could  decrease  the  maximum  error  by  a  factor  of  one-half  or  more, 
especially  when  the  comparisons  arc  to  be  made  between  a  number  of  data 
points  all  taken  on  the  same  channel  and  in  the  same  manner. 


APPKiUlX 


Data  £or  various  curves  are  co-rpiled  in  the  Appendix  and  are  identified 
with  the  particular  curve  by  figure  nurhcrs. 


Roc*  Velocity 
ft/spc 

1018  STEZL  DATA 
Flfi;-  1.15 

Strain 
i  /  £  n 

Rod  Velocity 
ft/sec 

St  lain 
in/in 

12 

-  00?-4?^ 

116 

.00224 

14 

.0004^5 

131 

.00316 

29 

.0'.;:  >80 

141 

.00275 

43.5 

.00135 

143 

.00275 

53.4  . - 

.00157 

145 

.00291 

63 

.001795 

160 

.  0024 

. 00245 

63.1 

.00176 

180 

.00254 

66 

.0018 

If'O 

.00246 

68 

.00198 

188.5 

.00245 

72 

.0019 

.00201 

i>0.5 

.0022 

75.5 

.0020 

191 

. 00242 

76 

.0024 

198 

.00251 

76.7 

.00238 

202 

. 00246 

79.5 

.00256 

203 

.00240 

83 .  ^ 

.00237 

205 

. 00243 

88 

.0024 

206 

.0025 

90 

.00251 

207 

,00246 

95 

.0027 

208 

,0024 

95.5 

.00275 

220 

.00232 

96 

.00254 

1 

223 

. 00234 

101,5 

.00223 

224 

.0022 

. 00264 

104 

.0023 

110 

.00287 

.00221 

TABLE  A-! 

1018  STEEL  DATA 
(cont  tn^ieti) 


Rod  Velocity 
ft /sec 

Stra in 
in/ir» 

237 

.0024 

247 

.00252 

248 

.00231 

.00259 

256 

.00243 

261 

.00253 

270 

.00237 

299 

.00237 

301 

.00252 

317 

.0030 

418  '  " 

.00251 

00251 


Rod  Velocity 

Srr.'T  •'r: 

2 

Rod  Velocity 

Strain 

ft /see 

in/  ir 

1 

5 1 /sec 

i  n  /  i  n 

34,5 

,  *  I-  j  i.  ^ 

i 

l‘f 

40.3 

.00! 19 

I'll, 

.00315 

48,3 

.00138 

197.2 

52.6 

197.5 

53.5 

.00175 

201,5 

72.5 

.00204 

221 

, 00433 

74 

.00218 

236 

. 00396 

87 

.00254 

249 

.00485 

92 

.00261 

269 

.00275 

96 

.0026 

281 

.00366 

100.8 

.0031 

285 

. 00422 
.00376 

116 

.00259 
. 002  78 

297 

.00321 
. 00407 

120.7 

,00318 

340 

.00467 

131 

.00314 

525 

.00415 

132 

.00324 

133 

.00315 

.00276 

138.5 

.00341 

141.5 

.00376 

158 

,00311 

1()0 

. 0044  7 

163 

.00442 

1  Velocity 
ft/sec 

TABLE  A- 3 
COPPER  DATA 

Fig.  1-19 

Strain 
in/  in 

Rod  Velocity 
ft^sec 

Strain 

in/in 

10.7 

.0051 

126.3 

.00298 

.00205 

13.2 

.00<K> 

138 

.00292 

21.8 

.000976 

. 00243 

28.8 

.001_6 

182 

.00268 

31.9 

.0013* 

191 

.00221 

.00321 

32.9 

.00135 

33.5 

.0017. 

33.7 

.001. ,5 

43.75 

.00147 

52.5 

.0;)19 

54.75 

. 00202 

66.2 

.00219 

73.2 

.00191 

.00209 

79.6 

.0023 

85.5 

.00231 

86 

.00231 

91 

. 00248 

91.8 

.00316 

92.5 

.00265 

103 

.00255 

113 

.00255 

122 

.00216 

124.6 

.00333 

80 
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